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Now  define  a  new  variable  *'*%  and  let  4*/’s*/;/rtr^then  the  integrals  can  be 
written  in  terms  of  dimensionless  quantities  as  follows: 

We  first  combine  the  elastic  and  plastic  cases  in  one  set  of  integrals 
noting  that 

a.  For  the  elastic  case  A-o 

b.  For  the  plastic  case  J  * 

The  combined  integrals  are  then 
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where 


(over) 


Thus  the  integrals  are  dimensionless  quantities  which  are  functions 
of  the  dimensionless  ratios 

Equation  [40]  is  the  form  that  the  electronic  computer  will  use 
to  compute  the  integrals  numerically  for  a  given  shape 
given  parameters  and  a  series  of  values  of 


It  should  be  noted  that  all  expressions  in  the  report  pertaining  to 
the  plastic  region  are  correct  if  <2^  .  This  condition  holds 

for  a  large  number  of  cases  in  view  of  the  fact  that  we  are  concerned 
with  large  deformations  and  large  plastic  strains.  If  this  condition 
is  relaxed  then  expression  [10]  must  be  rewritten  as 


Glcfei  -hj ^ Oi<s/e:J dZ. 
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The  last  term  must  be  subtracted  from  all  \/  expressions  if  we  are 
dealing  with  small  plastic  deformations  (  i.e.  if  £l;  «  Cj  ) 

In  the  example  worked  out  in  the  report,  eq.  [32]  would  then  become 


so  that  the  new  term  is  negligible  if 
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For  very  thin  shells  undergoing  very  large  plastic  deformations  if  X  s  / 
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J  J  ±-Jh  A  \  J«77d)  c*c/<fc/x-  (see  [16]) 
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For  some  physical  cases  *'%**!  %*<■  /  so  that 
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Thus  for  the  perfectly  plastic  case 
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ABSTRACT 


A  theory  for  the  dynamic  elasto-plastic  deformations  of  thin  shells 
is  presented.  The  deformation  is  obtained  by  equating  the  energy 
input  to  the  energy  absorbed  by  the  structure.  Large  deflections 
and  strain  hardening  are  considered.  The  theory  is  applied  to  the 
case  of  a  shell  undergoing  axially  symmetric  collapse. 
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I.  Introduction 

Shell  theory,  especially  the  dynamic  plasticity  of  shells,  is 
in  such  a  state  at  the  present  time  that  one  can  only  hope  to 
obtain  approximate  solutions  to  the  problems  that  are  of  prac¬ 
tical  interest  today.  Very  little  work  has  been  done  on  the 
dynamic  plasticity  of  shells  due  in  part  to  the  mathematical 
complexities  of  the  theory  and  probably  due  to  the  lack  of 
experimental  evidence  with  which  to  check  the  results  of  the 
theoretical  developments.  Fortunately  the  electronic  computer 
will  enable  us  to  overcome  some  of  the  mathematical  difficulties. 

Extensive  use  is  being  made  of  shells  in  missiles  and  submar¬ 
ines.  Therefore  more  complete  experiments  are  being  conducted 
on  these  structures  and  more  experimental  evidence  is  becoming 
available  to  those  working  in  the  field  of  dynamic  elasticity 
and  plasticity.  More  theoretical  development  on  this  problem 
is  needed  and  it  is  for  this  purpose  that  the  present  report 
has  been  written.  This  report  will  contain  an  approximate 
method  for  predicting  large  elastic  and  plastic  deformations 
of  shells  under  static  and  dynamic  loads. 

II.  Theory 

A.  Energy  relations 

Assume  that  the  shell  is  exposed  to  an  impulsive  load  of 
short  duration  which  imparts  an  initial  velocity  to  the 
structure.  The  problem  is  then  equivalent  ^o^one  in  which 
the  structure  has  an  initial  Kinetic  energy.  '  If  the 
energy  absorbed  by  the  shell  is  plotted  against  the  lateral 
deflection  in  both  the  elastic  ard  plastic  regions,  then 
according  to  a  previous  reference^  there  can  be  a  region  of 
instability  in  the  elastic  region  which  corresponds  to  snap 
buckling  of  the  structure.  If  no  instabilities  occur  in 


Witmer,  et  al  "Response  of  Plates  and  Shells  to  Intense 
External  Loads  of  Short  Duration, "  WADD  Technical  Report 
60-433,  Apr.  1960  (Mass.  In.  Tech.) 

“Tian,  T.  H.H.,  "Dynamic  Response  of  Thin  Shell  Structures," 
Plasticity,  Proc.  Sec.  Symp.  on  Naval  Struct.  Mech . ,  1960. 
p.  443. 

Radkowski,  P.P.  et  al,  "Studies  on  the  Dynamic  Response  of 
Shell  Structures  and  Materials  to  a  Pressure  Pulse, "  Avco 
Corp.,  1961,  AFSWC-TR-61-31  (II) 
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the  elastic  region,  when  the  structure  goes  plastic  there  could  be 
a  plastic  instability  point.  If  this  point  does  not  exist  and  the 
structur  is  loaded  further  then  there  should  be  a  point  at  which 
a  slight  increase  in  energy  will  correspond  to  a  relatively  large 
increase  in  deflection.  This  should  happen  when  the  structure  is 
about  to  collapse. 


The  initial  kinetic  energy  can  be  determined  in  terms  of  the  applied 
impulse.  This  impulse  will  in  general  turn  out  to  be  a  nonlinear 
function  of  the  lateral  displacement  if  large  strains  occur  and  if 
the  corresponding  large  deflection  expressions  are  employed.  It 
will  be  assumed  that  all  the  initial  kinetic  energy  that  is  imparted 
to  the  shell  by  the  impulse  goes  into  energy  of  deformation.  The 
final  deformation  will  occur  when  the  shell  comes  to  rest  so  that 
the  final  kinetic  energy  is  zero. 


Schuman  has  conducted  experiments  on  a  large  number  of  aluminum 
and  ;■»  eel  cylinders  subject  to  nonuniform  blast  loading.  He  found 
two  characteristic  types  of  failure  which  are  shown  in  Figure  1. 


!  B/ost 


Fig.  la.  Failure  Type  1 


Fig.  lb.  Failure  Type  2 


Unpublished  data  on  tests  performed  by  William  Schuman  at  the 
Ballistic  Research  Laboratories  at  Aberdeen  Proving  Ground. 
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It  seems  that  failure  type  1  is  a  collapse  with  a  failure  hinge  at 
the  center  of  the  cylinder.  Failure  type  2  resembles  a  buckling 
type  failure.  Whether  or  not  these  types  of  failure  fit  into  the 
catagories  of  plastic  collapse  and  buckling  instability  remains  to 
be  seen.  It  is  proposed  to  use  these  general  failure  shapes  found 
experimentally  as  the  deformation  patterns  in  the  energy  expression. 
In  order  to  determine  the  final  plastic  deformation  due  to  a  given 
impulse  the  energy-deflection  relation  will  first  be  computed.  The 
initial  kinetic  energy  will  then  be  equated  to  the  energy  of  de¬ 
formation  absorbed  by  the  shell.  The  applied  impulse  can  be  written 
in  terms  of  the  initial  kinetic  energy  as  will  be  shown  later  in 
the  report.  The  final  result  will  be  a  curve  of  applied  impulse 
as  a  function  of  final  deflection.  This  curve  will  implicitly  con¬ 
tain  all  the  nonlinearities  in  the  strain-deflection  relations  and 
all  the  details  of  the  elastic-plastic  stress-strain  law.  The  re¬ 
gions  of  buckling  and  collapse  can  be  determined  from  the  energy- 
deflection  curve  as  shown  in  Fig.  2 


Fig.  2  Buckling  and  Collapse  Curves 

The  deformation  energy  (or  work  done  by  the  internal  forces)  per 
unit  volume  of  an  elastic-plastic  body  can  be  written^ 

Iliouchine,  "Plasticite  (Translated  from  the  original  Russian  into 
French),  Edition  Eyrolles,  1956,  p.  98. 
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The  curve  of  »/3.  describes  the  stress-strain  law  of 

the  material  as  shown  in  Fig.  3 


Assuming  an  incompressible  material  ( &  -  o )  and  particularizing 
our  analysis  to  a  thin  shell 


The  notation  is  as  shown  in  Fig.  4 


Fig.  4  The  Shell 

Ihe  expressions  for  the  strains  are^ 

£ 

Love,  "Mathematical  Theory  of  Elasticity,"  Fourth  Edition,  Dover 
Publications,  1944,  p.  529. 
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-y-2ZT 


[3] 


where  €tj  )r  are  the  midsurface  strains,  j?  is  the  radial 

distance  from  the  midBurface  to  any  element  within  the  thickness 
of  the  shell  and  AT  A".  T  are  the  curvatures  and  twist.  The 

9  J 

values  of  ^  J"  are  (allowing  large  deflections)  ' 


1 
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/<2±f  \ r 

(  o%> J 


2  £k  ( &Jcf  *  C*- 
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and  the  curvatures  and  twist  are 


8 


hr  = 

I 


/fTr  Of  +±  JyT 

1  ax  icfv  a x 


[5] 


^  £Lit 

*  c)*. 


A  general  deformation  type  elastic-plastic  stress-strain  law  can  be 
written^  (assuming  a  Poissons  ratio  of  \) 


where 


C5-  f  %  <Vi 

W4. 


[6J 
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7 

Fung  and  Sechler,  "Instability  of  Thin  Elastic  Shells,"  Proc.  of  the 
First  Symp.  on  Naval  Struct.  Mech . ,  1958,  p.  118. 
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in  which  cue*i)=  o  for  the  elastic  region.  Now  neglect  all 
terms  in  the  strains  which  contain  c/  and  ,  assuming  that 

they  are  small  compared  to  ajlT  and  its  derivatives* 

Then 


_ /  /  Jttr  \ 1 

~  2.  (  OX.  J  “ 


Ojcx 


-  Quri+c 

Ox  ’ 
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j. 

* 


[8] 


Consider  an  elastic-linear  hardening  incompressible  material 
which  has  a  stress  strain  law  as  shown  in  Fig.  5 


♦This  assumption  seems  to  have  been  excellent  in  predicting  very 
large  strains  in  panels,  see  J.  Greenspon,  J.  Aero/Space  Sciences, 
27,  5,  May,  1960. 

♦♦Aluminum  and  steel  both  come  close  to  being  linear  hardening 
materials  -  see  ANC-5  "Strength  of  Metal  Aircraft  Elements,"  U.  S. 
Government  Printing  Office. 
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It  is  planned  to  use  this  stress-strain  law.  The  energy-deflection 
relation  will  then  be  computed  and  a  search  will  be  made  for  regions 
of  instability  in  the  elastic  and  plastic  regions  and  a  region  of 
collapse  in  the  plastic  regime.  If  the  value  of  energy  correspond¬ 
ing  to  instability  is  smaller  than  the  amount  corresponding  to  col¬ 
lapse,  the  structure  will  buckle  first?  however  if  the  collapse  energy 
is  smaller  than  the  energy  corresponding  to  instability,  the  struc¬ 
ture  will  collapse  first. 

Going  back  to  the  general  expression  for  the  energy,  and  substituting 
the  stress  strain  law  of  eq.  [9] 


[10] 


[11] 


Substituting  the  expression  for  <Zi  from  [2] 


[12] 
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Using  the  strain  expressions  [8]  and  integrating  with  respect  to  2 


f  c/< . 


[13] 


The  strain  terms  in  equation  [8]  which  involve  2  are  associated 
with  bending,  the  other  terms  are  the  membrane  strains.  For  very 
thin  shells  sometimes  the  bending  strains  can  be  neglected  and  the 
resulting  energy  expression  will  be* 


a  ctyc/x~ 


[15] 


Some  materials  can  be  considered  perfectly  plastic.  For  this  case 
X  =  /  in  the  plastic  region  and 


3'  ^ 


[16) 


♦Obtained  by  setting  -O  in  [13J;  note  that  L'Hopital’s 

rule  must  be  applied  to  the  second  bracket. 
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The  energy  of  deformation  in  the  elastic  region  can  be  obtained  from 
[13]  by  letting  X  =  O  •  However  a  more  accurate  expression  for  the 
elastic  region  can  be  obtained  by  starting  with  the  elastic  stress 
strain  law  for  any  Poissons  ratio,  i.e. 


^  £*■ ) 

x 4—  w 

if  y 


The  elastic  energy  is 


1 12a] 


Thus 


^  4r <?)]* dcp 

w  A  f 


[13a] 


where 
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a*-  2.  L  (  *<)0S 

/<r  '¥J  (  a*  j*'  3fi>  +  *+(  -a7‘(jZ3?J 


[14a] 
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In  some  cases  the  deformation  of  the  shell  may  be  axisymmetr ic  and 
the  deflection  -uf~  will  then  be  independent  of  .  For  this  case 

expression  [13]  becomes 


[18] 


This  concludes  the  basic  energy  equations.  The  relation  between 
the  kinetic  energy  input  and  the  applied  impulse  will  now  be  con¬ 
sidered. 


B  Impu] se-Energy  Relation 

Let  I  be  the  impulse  per  unit  mass  applied  to  the  shell.  The  impulse 
momentum  relation  for  an  elemental  mass  c!r*\  can  be  written 


[19] 


where  *uJr  is  the  lateral  velocity  imparted  to  the  mass  by  the  impulse* 
Thus 


-  X  [20] 


The  kinetic  energy  imparted  to  the  shell  is 


*Only  the  lateral  velocity  is  being  considered, 

being  neglected. 


u  and  'V  are 


[21] 


T=  f  £  -  f  fj>JU:W 

where  is  the  mass  per  unit  area  and  </*  is  an  elemental  area 

The  impulse  can  vary  over  the  surface  of  the  shell  in  the  same  way 
that  the  pressure  can  vary  over  the  surface.  Therefore  write  e 

impulse  as 


Tfx,$)  =  t0 


[22] 


Thus  from  [21] 


[23] 


Equating  the  initial  kinetic  energy  to  the  energy  of  deformation  ab 
sorbed  by  the  shell  the  expression  for  the  impulse  per  unit  mass  be 

comes 


2. _ 


[24] 


The  total  impulse  on  the  shell  will  then  be 
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It  should  be  emphasized  that  the  expressions  for  the  impulse  hold 
when  the  shell  has  reached  its  final  state  of  deformation  and  has 
come  to  rest  or  when  an  unstable  point  is  reached  where  the  kinetic 
energy  is  zero. 

C  Static  loading 

In  the  case  of  static  loading  the  energy  of  deformation  is  equated 
to  the  work  done  by  the  applied  load.  The  result  will  be  a  curve  of 
load  vs  deflection.  The  static  buckling  or  collapse  load  can  then 
be  obtained  in  a  similar  fashion  as  described  bdrore;  i.e.  first  the 
energy  deflection  relation  is  plotted  and  the  deflection  at  which 
instability  or  collapse  occurs  is  determined  as  shown  in  Fig.  2.  The 
load  corresponding  to  this  deflection  then  defines  the  static  buck¬ 
ling  or  collapse  load.  Since  the  energy-deflection  relation  will  in 
general  be  nonlinear  and  probably  will  not  contain  sharp  discontin¬ 
uities,  the  collapse  load  will  not  be  well  defined.  However  it  is 
expected  that  a  region  of  collapse  can  be  determined.  The  work  done 
1  j  the  static  load  will  be 


W  -  J  p  )  /LAS' 


[28] 


where 


is  the  applied  pressure  and  the  lateral  deflection. 
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D.  A  simple  application 

Although  extensive  calculations  will  be  done  in  the  near  future  using 
the  relations  derived  here,  a  very  simple  application  will  nott  be  con¬ 
sidered  in  which  the  integrals  can  be  evaluated  analytically.  Con¬ 
sider  the  axially  symmetric  lateral  loading  of  a  perfectly  plastic 
shell  subjected  to  both  static  and  dynamic  loading  which  is  assumed 
uniform  over  the  shell.  The  problem  is  to  estimate  the  static  col¬ 
lapse  load  and  to  determine  the  impulse-deflection  relationship  in 
the  plastic  region. 

It  has  been  found  experimentally^  that  cylindrical  shells  under 
axially  symmetric  external  load  usually  collapse  with  a  deformation 
pattern  resembling  two  frustrums  of  a  cone  with  a  hinge  circle  at 
the  center  of  the  cylinder  as  shown  in  Fig.  6 


Fig.  6  Collapse  Pattern  of  Cylindrical  Shell 
under  Axially  Symmetric  Load 

According  to  the  sign  convention  used  in  this  analysis  aaJT  is  positive 
inward.  The  axes  and  deformation  pattern  are  shown  in  Figure  7. 


Fig.  7  Deformation  Pattern  of  Shell 


Thus 


Integrating 


[291 


4 


[30] 


Assume  that  the  shell  and  deformation  pattern  are  such  that* 


a 


<  c  / 


>  >  / 


Thus  the  energy  expression  reduces  to 


[31] 


[32] 


For  static  loading  the  work  done  by  the  external  uniform  lateral 
x oct u  is 


"This  will  be  true  of  the  majority  of  practical  cases. 
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[33] 


Z - 


Equating  the  internal  energy  absorbed  by  the  shell  to  the  work  done 
by  the  external  load  it  is  found  that 

'P<‘~  ar£.  [34] 

Expression  [34]  implies  that  the  load  will  not  be  dependent  on  the 
deflecti throughout  the  plastic  region.  Therefore  this  load  will 
correspond  to  any  deflection  in  the  plastic  region  and  can  therefore 
be  termed  the  static  collapse  load. 

It  is  interesting  to  note  that  the  load  corresponding  to  yield  in  an 
elastic  shell  with  ordinary  hoop  tension  is 

Gfi)t  -  °IJ:  [35] 


Thus  [34]  predicts  that  the  static  collapse  load  is  about  15%  higher 
than  the  load  at  which  yielding  will  start. 

For  shells  in  which  is  not  much  greater  than  unity  higher  order 

terms  in  the  energy  must  be  included.  If  second  order  terms  in 
are  retained  in  [30]  the  energy  becomes 


vf  (  -zz:  +  +  / 


[36] 


Therefore  the  static  load-deflection  relation  in  the  plastic  region 
will  be 


i>.  =  (•  *  j  ^4-  ) 

/J  ^  l  "  'i.  "v*.  J 


[37] 


This  will  give  a  load  deflection  curve  which  looks  schematically  as 

wn  in  Fig.  8. 
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Fig.  8  Load  Deflection  Curve  in  Plastic  Region 
Including  Nonlinear  Terms 

The  effect  of  the  nonlinear  terms  is  to  stiffen  the  shell  so  that 
the  plastic  deflection  for  a  given  load  will  be  a  definite  value. 
For  the  linear  case  the  collapse  load  would  produce  any  plastic 
deflection. 

Using  relations  [32]  or  [36]  with  [27]  the  impulse-deflection  curve 
in  the  plastic  region  can  be  obtained.  A  schematic  of  the  results 
using  these  two  energy  relationships  is  given  in  Figure  9. 
with  [32] 


(ZoMf-  mijijeA. 

v/T  7To£ 


[38] 


with  [36] 


]/3  i  2-c*.  +  J  m*jL 


139] 
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Fig.  9  Impulse  Deflection  Relation 

The  static  collapse  load  does  not  define  collapse  in  a  dynamic  prob¬ 
lem.  In  fact  if  the  load  is  applied  dynamically  the  actual  magni¬ 
tude  of  the  load  can  be  much  greater  than  the  static  collapse  load 
and  still  not  result  in  plastic  deformation.  It  is  the  load-time 
relation  that  is  important  and  the  impulse  is  a  lumped  parameter 
which  essentially  is  a  measure  of  the  load  time  history. 

III.  Discussion 

The  great  advantage  of  the  energy  absorption  method  as  described 
here  and  elsewhere  in  the  literature1'  is  that  many  complications 
can  be  considered  in  the  analysis  and  still  result  in  a  tractable 
problem.  For  example  it  is  seen  that  we  have  considered  here  non¬ 
linearity  in  the  strains,  strain  hardening  in  the  stress-strain 
law  and  nonuniform  dynamic  loading.  This  type  of  theory  has  one 
main  disadvantage  -  we  have  to  assume  a  deflection  pattern.  How¬ 
ever  it  must  be  realized  that  even  in  the  more  sophisticated  plas¬ 
ticity  theories  assumptions  on  the  pattern  must  also  be  made. 

Based  on  careful  study  of  the  existing  plasticity  flow  theories  of 
today  it  is  difficult  for  this  author  to  see  how  the  important  com¬ 
plications  such  as  nonlinearity  in  the  strains  and  stress-strain  law 
in  addition  to  nonuniform  dynamic  loading  can  be  included  and  result 
in  a  solvable  problem.  In  such  plastic  dynamic  biaxial  stress  prob¬ 
lems  this  author  sees  very  little  hope  for  any  but  the  deformation 
theories  and  the  energy  type  procedure  for  actually  getting  reasonable 
answers  in  a  reasonable  time. 
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